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ABSTRACT 

For measurement of the AGN luminosity function and its evolution, X-ray selection samples all 
types of AGN and provides reduced obscuration bias in comparison with UV-excess or optical surveys. 
The apparent decline in optically-selected quasars above z '^ 3 may be strongly affected by such a 
bias. The Chandra Multiwavelength Project (ChaMP) is characterizing serendipitously detected X- 
ray sources in a large number of fields with archival Chandra imaging. We present a preliminary 
measure of the co- moving space density using a sample of 311 AGN found in 23 ChaMP fields (~ 1.8 
deg^) supplemented with 57 X-ray bright AGN from the CDF-N and CDF-S. Within our X-ray flux 
(/o. 3-8.0 kev > 4 X 10"^^ erg cm~^ s~^) and optical magnitude (r' < 22.5) limits, our sample includes 
14 broad emission line AGN at z > 3. Using this X-ray selected sample, we detect a turnover in the 
co-moving space density of luminous type 1 AGN (log Lx > 44.5; units erg s""'^; measured in the 0.3- 
8.0 keV band and corrected for Galactic absorption) at z > 2.5. Our X-ray sample is the first to show 
a behavior similar to the well established evolution of the optical quasar luminosity function. A larger 
sample of high redshift AGN and with a greater fraction of identified sources, either spectroscopic or 
photometric, at faint optical magnitudes (r' > 22.5) are required to remove the remaining uncertainty 
in our measure of the X-ray luminosity function, particularly given the possibility that AGNs might 
be more easily obscured optically at high redshift. We confirm that for z < 1, lower luminosity 
AGN (log Lx < 44.5) are more prevalent by more than an order of magnitude than those with high 
luminosity. We have combined the Chandra sample with AGN from the ROSAT surveys to present 
a measure of the space density of luminous type 1 AGN in the soft X-ray band (0.5-2.0 keV) which 
confirms the broad band turnover described above. 
Subject headings: galaxies: active — quasars: general — X-rays: galaxies — surveys 



p ,■ 1. INTRODUCTION 

Q . Optical sur veys h ave measured the evolution of QSOs 

^ out to 2 ~ 6 l|Fan e t al. 2004). The most dramatic fca- 

jy-j ture found is the rise and fall of the co-moving space 

Ci ' density with peak activity at z '^ 2.5. A systematic 

L^ ' decrease in luminosity from z '^ 2 to the present is evi- 

. ,_( ' dent with very few intrinsically bright QSOs in the local 

^^ universe (e.g.. Groom ct al. 2004). This fading of the 

Vh QSO population is attribut ed to a decre ased fuel sup- 

p^ , ply and/or fueling ra te fe.g. .jCavaliere fc Vittorini _2000: 

iKauffmann fc Haehn clt 20 Q3) . The dropoff in the space 



density at z > 3 (e .g., Warren. He wett fc Osmer I 



1^ 



3tt — 

[ Schmidt. Schneider & Gunn 1995; iFan et al. F I200 
|Wolf ct al. 2003) could represent the growth phase of su- 
permassive black holes (SMB Hs). possibly under a veil 
of obscuration ( Fabian 11999(1 . 

It has become clear in the past decade that X-ray selec- 
tion of AGN offers many benefits over detection methods 
in other wavebands. X-rays can penetrate large absorb- 
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ing columns of gas which can effectively hide an accreting 
black hole in the optical. The existence of a significant, 
missed population of obscured AGN is required to ex- 
plain the spectral shape of the Cosmic X-ray Background 
(CXRB; e.g.,rGiUi. Salvati & Hasinger 2001). The cur- 
rent generation of X-ray observatories ( Chandra, XMM- 
Newton) are probing the faint source p opulation respon- 
sible for the bulk of the CXRB (e.g.. [Alexander et al. I 
l2003lm.osa,ti et al. 1120021 iHTsinger et al. 112001: 

X-ray surveys with Chandra, XMM-Newton and 
ASCA are refining our knowledge of the AGN X- 
ray luminosity function. The peak, in the co-moving 
space density of low luminosity (i2.o-8.o kcV < 10*"' 
erg s~^) AGN found i n hard (2-10 keV) X-ray surveys 
JBarger et al. I l200| iCowie et al. I l200a iFiore et aL I 



T lSteffan et al. 20031: lUeda et al. II2003D . occurs at 
z ~ 1. This contrasts starkly with the behavior of 
the luminous (L2.0-8.0 kcV > 10** erg s~^) QSOs which 
are most prevalent at z ~ 2.5. Using a highly com- 
plete sample of 941 AGN selected in the soft band, 
iHasinger. Mivaii fc Schmidt! l|2005|) show the same lu- 
minosity dependence and extend the space density mea- 
surement of low luminosity AGN {Lx < 10** erg 
s~*) out to z ~ 3. These latest results clearly re- 
quire a luminosity dependent density evolution model 
l)Mivai i. Hasinger fc Sc hmidt I 12000^ in contrast to a 
'pure' luminosity evo lution model, descr iptive of opti- 
cally selected QSOs llCroom et al. II2004D. With these 
results, semi-analvtic models (jMenci et al71 120041) posit 
that the more luminous QSOs illuminate early epochs 
when most of the massive galaxies are forming (z > 2), 



Silverman et al. 



thereby inducing high accretion rates, whereas lower lu- 
minosity AGN dominate at a later period (z < 2) when 
most of the galaxies have fully assembled. 

X-ray selected QSOs from ROSAT have hinted at 
a constant space dens ity between 1.5 < z < 4.5 
IjMivaii. Hasinger fc Schmidt 20CKI ) for the most lumi- 
nous AGN (log Lg. 3-8.0 kev > 44.5). The lack of a de- 
cline in the space density at z > 3, in contrast to the 
behavior seen in optical surveys, could be evidence for 
obscuration at early epochs. We can test this model 
(Fabian 1999) by measuring the luminosity function and 
co-moving space density of X-ray-selected AGN over a 
wide area to compile a significant sample at z > 3. Un- 
fortunately, the ROSAT sample includes only 7 QSOs 
at these high redshifts. Of these, only 5 QSOs have 
luminosities high enough to be detected over a broad 
range of redshift (0 < z < 5) and thereby included 
in the measurement o f the co-m oving space density 
IjMivai i. Hasinger & Sc hmidt 1 120001 . While the Chan- 
dra and XMM-Newton Deep field observations have great 
scientific merit, their narrow area provides only 4 high 
luminosity (log Lq. 3-8.0 koV > 44.5) AGN at z > 3. Re- 
cent measurernent of the hard X-ray lum i nositv function 



cm ^ s ^), we effectively smooth out any effects from 
large scale structures such as those evident in the CDF- 



la rge scale stru ctures 
S llGilli et al. II2OOI 



llBarger et aTl I200I iFiore et al. 1 12001 iSteffan et all 
120031 lUeda et al. Il2003|) have limited numbers of AGN 
at these high redshifts. A wider area survey is needed 
to compile a significant sample of highly luminous QSOs 
with redshifts greater than 3. 

We present a preliminary measure of the co-moving 
space density of X-ray selected AGN out to z ^ 4 
from the Chandra Multiwavelength Project (ChaMP). 
We supplement the ChaMP sample with AGN of com- 
parable X-ray and optical fluxes from the Chandra Deep 
Field North (CDF-N) and South (CDF-S) to boost our 
spectroscopic completeness at faint optical magnitudes 
and include one additional z > 3 QSO. We merge the 
subsample of these AGN detected in the soft (0.5-2.0 
keV) band with those from ROSAT to directly compare 
with Miyaji et al. (2000, 2001). A full presentation of the 
ChaMP X-ray luminosity function is forthcoming (Silver- 
man et al. 2005, in preparation). We assume Ho — 70 km 
s~^ Mpc~^, r^A = 0.7, and fiivi — 0.3 with the exception 
of Section O 

2. CHANDRA MULTIWAVELENGTH PROJECT 

The ChaMP IjKim et al. 112004311076^ et al. II200I is 
carrying out an extragalactic X-ray survey encompass- 
ing 10 deg^ using serendipitous detections in archival 
Chandra fields. Chandra's small point spread func- 
tion (~ 1" resolution on-axis) and low background al- 
low sources to be detected to fainter flux levels (^ 
10~^^ erg cm~^ s~^) than any X-ray observatory past 
or present, th us enabling the detection of high redshift 
(z ^ 5) AGN JTreister et al. II200I IBarger et al. MM 
iCastander et al. ll2003HSilverman et al. 1 12002^. For the 
present study, we use the full Chandra energy range (0.3- 
8.0 keV) to detect the absorbed sources missed by pre- 
vious optical, UV or soft X-ray surveys and take ad- 
vantage of the high collecting area at soft energies to 
detect the faint, high redshift AGN. At z > 3, we are 
sensitive to AGN with absorbing columns up to ~ 10^^ 
cm~^ due to the favorable fc-correction. With a sample 
of AGN selected from a large number of non-contiguous 



Chandra fields that reach similar depths 
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We have chosen 23 Chandra fields (1.8 deg^) for which 
we have acquired extensive foUowup optical imaging and 
spectroscopy. The deepest observations have exposure 
times that are sensitive to sources with /o.3-8.0keV > 
8 X lO"^*" erg cm"^ s~^ A full description of the ChaMP 
image reduction an d analy sis pipeline XPIPE can be 
found in 'K im et al."l l)2004b(l . With our 4m MOSAIC op- 
tical imaging, we are able to identif y counterparts to the 
Chandra sources down to r' ~ 25 l)Green et al. 1 12004'). 
We acquire optical imaging in three (q' , r^ and i') Sloan 
Digital Sky Survey (SDSS) filters (Fu kugita et al. 11 996). 
Optical colors provide preliminary source classification 
and crude photometric redshifts. These diagnostics are 
required to characterize the optically faint X-ray sources 
that cannot be identified with spectroscopy and in par- 
ticular, their influence on the X-ray luminosity function. 
The use of the SDSS photometric system allows more 
direct compar ison between the ChaM P and the SDSS 
AGN surveys IjSilverman et al. Il2004|) . Optical spectro- 
scopic followup currently focuses on identifying counter- 
parts with r' < 22.5 for which spectra can be acquired 
on a 4-6m (i.e. MMT, Magellan, WIYN, CTIO Blanco) 
class telescope. To date, we have spectroscopically clas- 
sified a sample of 358 AGN detected in the broad band 
(0.3-8.0 keV) in these 23 ChaMP fields. 

3. X-RAY SENSITIVITY AND AREA COVERAGE 

There are complications with measuring the X-ray lu- 
minosity function that must be put into context. The 
difficulty in focusing X-rays results in a point source func- 
tion (PSF) and flux sensitivity that varies across the field 
of view. The PSF degrades as a function of off-axis an- 
gle, decreasing the flux sensitivity. In terms of measuring 
the luminosity function, it is not a trivial task to deter- 
mine the incompleteness at the faintest X-ray fluxes and 
the actual sky area over which a source of a specific flux 
would be detected. In addition, the exposure times of 
these Chandra flelds range from 17 to 114 ksec generat- 
ing a wide range of limiting fluxes even on-axis. 

To characterize the sensitivity, completeness, and sky 
area coverage as a function of X-ray flux, a series of simu- 
lations were performed. The full details will be presented 
in an upcoming ChaMP X-ray analysis paper (M. Kim, 
D.-W Kim et al., in preparation). The simulations con- 
sists of three parts, 1) generating artificial X-ray sources 
with MARX (MARX Technical Manual^) and adding 
them to real X-ray images, 2) detecting these artificial 
sources by wavdetect and extracting source properties 
with XPIPE identically as performed for actual sources, 
and 3) estimating the sky area coverage by comparing 
the input and detected source properties as a function 
of flux. The simulations are restricted to specific CCDs 
for ACIS-I (10, II, 12, and 13) and ACIS-S (12, 13, S2, 
and S3). These CCDs are closest to the aimpoint for 
each observation. Sources far off-axis (O > 12') are ex- 
cluded since the flux sensitivity is low and the PSF is 
degraded. These simulations allow determination of cor- 
rections for the source detection incompleteness at faint 
flux levels quantifled in the first ChaMP X-ray analy- 

* http://space.mit.edu/CXC/MARX 
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Fig. 1. — Area coverage as a function of broad (0.3-8.0 keV; 
top) and soft (0.5-2.0 keV; bottom) band X-ray flux for 23 ChaMP 
fields. A sinootlied curve is generated since small scale variations 
or sharp features are a result of the limited number of simulated 
sources and not real. 



sis paper IjKim et al. Il2004b|) . In Figure ^ we show the 
sky coverage determined from the simulations using the 
broad (0.3-8.0 keV) and soft (0.5-2.0 keV) band source 
detections for 23 ChaMP fields. We are surveying a sky 
area of 1.8 deg^ for the brightest sources. The sky area 
falls below 0.1 deg^ at the faintest flux levels. 

4. AGN SELECTION 

In these ChaMP fields, we find a diversity of objects 
(AGN, clusters, galaxies, and stars), although 85% of 
them are attribu ted to an AGN ( Green ct al. 2004; 
[Silverman et al. 1 ^004) . We show the optical magni- 
tude (r') as a function of 0.3-8.0 keV X-ray flux for 
the sources detected in 23 ChaMP fields (Figure |2Jl. We 
only include sources with greater than 9.5 net counts. 
Although sources with counts as low as 2 may be sig- 
nificant detections as a result of Chandrda low back- 
ground, we use a higher count limit to ensure that we 
have well-measured X-ray fluxes within our sample. In 
Figure 121 we label sources as classified from our optical 
spectroscopy. Extragalactic objects with strong emis- 
sion lines (Wa > 5 A) are labelled as either Broad Line 
AGN (BLAGN; FWHM > 1000 km s"!) or Narrow Emis- 
sion Line Galaxy (NELG; FWHM < 1000 km s"!). The 
BLAGN are equivalent to type 1 AGN. We loosely refer 
to the NELG with Lx > lO''^ erg s'^ as type 2 AGN 
since optical spectra may not cover all the emission lines 
needed to confirm photoionization from a non-thermal 
source. Extragalactic counter-parts with weak emission 
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Fig. 2.— Optical magnitude (r') vs. X-ray flux (0.3-8 keV) 
of X-ray sources found in 23 ChaMP fields. Optical spectroscopic 
classifications are indicated with the sample size. X-ray sources 
with no optical counterparts are shown by an arrow placed at the 
magnitude for a 5a detection from our optical imaging. 



line (Wa < 5 A) or pure absorption line spectra are clas- 
sified as Absorption Line Galaxy (ALG). In addition, a 
handful of stars are identified at faint X-ray fluxes. Six 
clusters have been found based primarily on their ex- 
tended X-ray emission. 

We calculate the rest frame (0.3-8.0 keV; erg s~^) X- 
ray luminosity for each extragalactic source, after cor- 
rection for Galactic absorption. The conversion from 
X-ray count rate to flux units (erg cm~^ s~^ ) is deter- 
mined from simulated detections on each CCD of a source 
with a powerlaw spectr um^ (/b cic ii^~ ^'"~^^;r = 1.7) and 
galactic absorption (Dickev fc Lockman 1990) . The ob- 
served luminosity is converted to the rest frame assuming 
a powerlaw spectrum with the photon index (F) set to 
the average value for that object type (BLAGN, NELG, 
ALG) based upon our X-ray spectral fit results (Aldcroft 
et al. in preparation). We find that the BLAGN have a 
mean spectral index (< F >) of 1.9 while the NELG have 
< F >=0.9 and ALG have flatter slopes < F >=1.5. We 
do not apply a correction for intrinsic absorption due 
to the uncertainty in the absorbing columns for most of 
these AGN as a result of their low X-ray counts. 

To construct a pure AGN sample, we require the de- 
rived (rest frame 0.3-8.0 keV) luminosity to exceed 10'*^ 
erg s~^ thereby excluding any sources that contain a sig- 
nificant stellar or hot ISM component. As shown in Fig- 
ure 121 we see that 92% of the extragalactic sources in 
the ChaMP not associated with an extended X-ray emit- 
ting cluster satisfy this criterion. We currently have not 
identified any NELG and ALG above a redshift of 1 in 
the ChaMP. Our optical spectroscopic limit (r' ^ 22) 
precludes classification of objects with heavy optical ex- 
tinction beyond this redshift, since a ~ 5L* elliptical 
host galaxy has r' ~ 23 at z ~ 1. We detect many 
optical type 2 AGN with Lx < lO'*'^ erg s'^ but no lu- 

et al. ll2004bl) provides energy con- 
version factors (EOF) for two models with F = 1.7 and F = 1.4. 



The ChaMP XPIPE 
3ion factors (EOF) for 
We chose the former since the photon index more closely resem- 
bles the majority of the X-ray source detections with spectroscopic 
identification. 
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objects found in 23 ChaMP fields (with the exception of 6 clusters). 
The dashed line marks the minimum luminosity required for the 
AGN sample. Optical spectroscopic classifications are shown as 
either BLAGN, NELG or ALG. 



minous (Lx > 10'''' erg s^-') optical type 2 QSOs. We 
do find 6% of the BLAGN to be X-ray absorbed with a 
substantial fraction having reddened optical colors which 
may prevent them f rom being selected in optical surveys 
IjSilverman et al. 120 04). These red BLAGN may be sim- 
ilar to those found by the Two Micron All Sky Survey 
(^rancis. Nelson & Cutri 2004) that ten d to have sig- 
nifica nt X-ray absorbing columns of gas IjWilkes et al. I 
|2002() . Deep optical foUowup of a broadband X-ray- 
selected sample will help determine the fraction of such 
red quasars to the space density of X-ray luminous AGN. 

5. SURVEY COMPLETENESS 

We have supplemented our C haMP sample with 
X-ray source s from the CDF-N (iBrandt et alJ 120011 
IBarger et al.l l200l and CDF-S ilSzokolv et al. I I200I 
with X-ray fluxes above our chosen limit (see below; 
/o. 3-8.0 kev > 4xl0~^^ ergcm"^ s~^). These samples al- 
low us to boost the number of spectroscopically identified 
sources at faint optical magnitudes (r' > 21) and include 
any additional z > 3 AGN. These are currently the only 
surveys with published optical spectroscopic identifica- 
tions at moderate completeness levels for sources with 
flux levels comparable to the ChaMP. In Figure 0] top, 
we plot the flux distribution of X-ray sources found in 
these 25 Chandra fields. The CDF-N and CDF-S clearly 
contribute most of the faint sources (log fx < —15). 
In Figure 01 bottom, we show the optical magnitude dis- 
tribution of counterparts to sources with X-ray fluxes 
above our chosen limit and those with optical spectro- 
scopic identifications. Optical magnitudes for the CDF 
N-l-S sources were converted from the Johnson (B, R) 
to the SDSS photometr i c syst em using the transforma- 
tion in lFukugJta et al.l l)199(i(l . By including these two 
deep Chandra fields, our area coverage is larger than that 
shown in Figure^ We add a fixed area of 0.22 deg^ since 
the area coverage in these two fields does not change 
above this bright flux limit by more than 3%. 

To evaluate the status of our optical foUowup and 
determine effective survey limits, we measure the frac- 
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Fig. 4.— top X-ray (0.3-8.0 keV) flux distribution of 1989 sources 
in the ChaMP and CDF N+S surveys (solid histogram). Those 
sources from only the CDF N-l-S are shown by the long dashed 
line. Those with optical spectroscopic identifications are shown by 
the short dashed line. The vertical line marks the X-ray flux limit 
of /o. 3-8.0 koV > 4 X 10~^^ erg cm~^ s~^. bottom Optical magni- 
tude (r') distribution of counterparts to X-ray sources with fluxes 
above our limit. Counterparts with spectroscopic identification are 
shown with a short (ChaMP) or long (CDF N-l-S) dashed line. The 
vertical line marks our chosen optical magnitude limit. 



tion of X-ray sources identified with optical spectro- 
scopic observations as a function of X-ray flux and op- 
tical magnitude. Since the X-ray sources do not fully 
sample the fx — t' plane (Figure 12)), we have imple- 
mented an ad aptive binning scheme identical to that 
of Sanders & F abian I (J2001) to generate a completeness 
map as a function of X-ray and optical flux for the com- 
bined sample (ChaMP -f CDF-N -f CDF-S). We create 
an array with 64x64 elements covering the /x — r' plane 
(Figure EJ ■ The identified fraction /id is measured at 
each element as the number of identified objects divided 
by the number of X-ray sources. If 10 or more X-ray 
sources are in a single element, then we assign /id as the 
completeness level. We then increase the bin size by a 
factor of two. If the binned element has 10 or more ob- 
jects, the individual elements not previously assigned an 
identified fraction arc set to this value. This procedure 
is iterated until the final binned element equals the size 
of the full array. In Figure [SJ we plot the final array (i.e. 
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Fig. 5. — Adaptively binned spectroscopic completeness map. 
The shading scale ranges from zero (lightest) to 100% (darkest) 
spectroscopically classified. Only spectroscopically identified AGN 
are shown with a black dot, those at 2 > 3 circled. Dashed lines 
show our chosen flux limits. A solid contour line has been over- 
plotted to mark the 30% completeness level. 



completeness map) as a grey scale image. 

Based on Figure we set the flux limits of our 
AGN sample. We choose a single X-ray flux limit at 
/o. 3-8.0 kcv > 4 X 10~^^ erg cm~^ s~^ that is brighter 
than the faintest X-ray detections in the ChaMP. This 
limit ensures that we minimize the bias towards optically 
bright AGN at the faintest X-ray flux limits. In order to 
limit our completeness correction to < 3.3, we set an op- 
tical magnitude limit at a level where at least 30% of 
the sources are identified, this leads to a magnitude limit 
of 22.5. With the remaining sample, the incompleteness 
level is used to correct the measurement of the space 
density (Section EJ. 

In Figure \^ we plot the location of AGN, includ- 
ing those from the Chandra Deep Fields in the /x — r' 
plane. Since a significant number of AGN fall in regions 
of low completeness at faint fluxes, we must exclude 8 
of the 22 AGN with z > 3 and log fx > -14.6 (cir- 
cled points). In particular, the high est redshift AGN 
(CXO MP J213945. 0-234655; z = 4.93: ISilverman et al. I 
I2002J) found by the ChaMP has been excluded not be- 
cause of its X-ray flux (/o. 3-8.0 kcV = 4.32 x 10^^^ 
erg cm~^ s""'^) but rather due to its optical magnitude 
(r' = 22.7) which falls below our limit. We do have a sig- 
nificant sample of 368 AGN above our flux limits with 14 
(13 from ChaMP) at z > 3, almost a factor of three larger 
than the high redshift sample available from ROSAT. In 
Figure we plot the Lx — z distribution of the total 
sample of 368 AGN selected in the 0.3-8.0 keV band. 

6. CO-MOVING SPACE DENSITY 

We implement the 1/Va method IjSchmidt II1968|) to 
generate an estimate to the co-moving space density (n 
in units of Mpc"'^; Equation^ in fixed redshift and lu- 
minosity intervals. For each L — z bin, Equation ^ is 
summed over all AGN (N) from our sample (Figure EJ) 
that fall within this bin. 
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As mentioned in Section O our survey has non- 
contiguous sky coverage with varying flux limits. To ac- 
count for this, we determine the accessible volume Va 
(Equation 121) over which an individual AGN i will be in- 
cluded in our sample, given the X-ray and optical flux 
limits. 
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dz 



dzdfl 
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We calculate the X-ray and optical limiting redshifts 



(^, 



7° 



) for each object. The integrand in Equa- 



tion |2| is summed from zi to Z2 where Z2 is the smaller 
°f ^lim ^^'i ^Hm if ^^^ AGN cannot be detected over 
the full redshift interval. The solid angle (dfi) is a 
function of X-ray flux as shown in Figure ^ plus an 
additional 0.22 deg^ contributed by the Chandra Deep 
Fields observations. A correction factor (/id), deter- 
mined from the location of an individual AGN in the 
fx — f' plane (Figure [SJ, compensates for the incom- 
pleteness in our spectroscopic identifications as detailed 
in the previous section. Each AGN contributes (/idKj)~^ 
to the space density in a specific luminosity and redshift 
interval. We estimate la errors based on a Poisson dis- 
tribution due to the small number of objects per red- 
shift bin. While there are known biases inherent in the 
1/Vn. method (e . g. M ivaii, Hasineer & Schmidt 200^ 
IPaee fc Carrera 1 120001) . the overall evolutionary trends 
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can be discerned (e.g. IBarger et al. 1120031 fCowie et al. I 
120031: iFiore et al. 112003(1 . A detailed X-ray luminosity 
function, will be presented in a future publication (Sil- 
verman et al. in preparation), while implementing more 
sophisticated analysis techniques (e.g. maximum likeli- 
hood method). 

6.1. Broadband (0.3-8.0 keV) 

We measure the co-moving space density (Equation^l 
of AGN detected in the broad band from 25 Chandra 
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Fig. 7. — Co-moving space density of AGN selected in the broad 
(0.3-8.0 keV) band. Vertical bars are estimates of the la error. 
Horizontal bars, centered on each data point, mark the redshift 
bin size. Horizontal bars with a downward arrow show the highly 
conservative upper limits that take into account the optically faint 
(r' > 22.5), unidentified X-ray sources. The upper limits corre- 
spond to the data which have the same line type connecting the 
points. 



fields. In Figure the number density is plotted in 
three luminosity bins (42.0 < log Lx < 43.0; 43.0 < 
log Lx < 44.5; 44.5 < log Lx < 46.0; L is in units of 
erg s~^) as a function of redshift. We must keep in mind 
that we are limited to the detection of obscured AGN 
(type 2) at luminosities below 10'^'* erg s~^ and redshifts 
below 1 due to our optical spectroscopic limit as shown 
in Figure |2I As a result, the highest luminosity range 
(44.5 < log Lx < 46.0) is restricted to type 1 AGN only. 
To account for the optically faint (r' > 22.5), spectro- 
scopically unidentified X-ray sources, we re move the opti- 
cal flux limit and implement the method of lCowie et al. I 
((2003} and Barger et al. (2003) to measure upper limits 
to the space density. We have assigned all unidentified 
sources to each feasible L — z bin with the redshift and 
luminosity set to the mean value. With this method, a 
specific source usually falls within multiple L — z bins. 

At z < 1, the low luminosity (log Lx < 44.5) AGN are 
more than an order of magnitude more numerous than 
those with high luminosity, thus confirming recent re- 
suits (ICowie et al. ll2003HFiore et a,1. ll20MlTTeda et al. I 
I2003D . This is evident even when considering the uniden- 
tified population since the upper limit to the space den- 
sity of luminous AGN (log Lx > 44.5) at z < 1 is fairly 
constrained. From Figure O we see that most of the 
sources contributing to the space density at z < 1 are 
non-BLAGN. These results agree with those reported by 
Eteffan et al. (2003) that show type 2 AGN as the dom- 
inant population at 0.5 < z < 1. The peak values of the 
two lower luminosity curves (n ~ 3 — 7 x 10~^ Mpc~^) 
agree with those from the soft (0.5-2.0 keV) X-ray lu- 
minosity function IjHasinger. Mivaii & Schmidt 120 05V ^" 
We find the space density a factor of ~3 less, for AGN 
with 43.0 < log Lx < 44.5 at z ^ 0.8, thaii the ha rd (2-10 
keV) X-ray survey results of lUeda et aTl (|2003() . These 

^^ Also reproduced in Fig ure 8a of a review article o n deep ex- 
tragalactic X-ray surveys by IBrandt fc Hasinger 1 120041) . 



comparisons seem reasonable since our broad band (0.3- 
8.0 keV) selected sample will be affected by absorption 
as known to occur in samples selected in the soft band 
(0.5-2.0 keV). 

For the luminous AGN (log Lx > 44.5), we measure a 
rise and fall of the number density with a peak at z ^ 2.5. 
With 14 AGN in the high redshift bin (3 < z < 4.2), we 
have evidence for a drop in the density that is inconsis- 
tent (at 3.4(t) with the peak value found at 2 < z < 3. 
This result is not strongly affected by our bright optical 
magnitude limit since most AGN at z > 3 do not all fall 
near the flux limit as shown in Figure O However, given 
the large numbers of unidentifled soures, the upper limits 
to the space density at z > 2 (Figure [Tj) show consider- 
able uncertainty. A higher fraction of identifled X-ray 
sources is required to determine if a decline is present 
with high signiflcance. The peak in the space density ap- 
pears to shift to lower redshift with decreasing luminosity 
as modelled by luminosity-dependent density evolutio n 
lITTeda et a,1. \mM iMhT aii. Hasi nger fc Schmidt \\2ml^ . 
in contrast to 'pure' luminosity evolutio n reported by 
the o ptical surveys of bright quasars (e.g. iCroom et al. I 
I200I . 

6.2. Soft band (0.5-2.0 keV): Chandra + ROSAT 

With Chandra's broad energy range, we can se- 
lect an AGN sample in the soft band (0.5-2.0 keV) 
to directly compare with the ROSAT results. The 
ChaMP AGN, supplemented by those in the Chandra 
Deep Fields, c an be co mbined with the ROSAT sam- 
ple (Mivaii. Hasinger fc Schm i dt 2000) to measure the 
number density of AGN in the soft band with signif- 
icant numbers of type 1 AGN up to z ~ 4. The 
ROSAT sample is a compilation of AGN found in sur- 
veys spanning the pencil beam (0.3 deg^), deep Lockman 
Hole ( Lchmann e t al. 2001) to the wid e area (2.0 x 10^ 
deg^) , shallow ROSAT Bright survey (jSchwope et al. I 
|2Q03). To directly compare with the published re- 
sults fro m the ROSAT and optical surveys such a s 
the 2dF (Groom et al. "20041) , SPS S (jFan et al. Il200l|) . 
and COMBO-17 (Wolf ct all 1200^ . we use cosmologi- 
cal parameters flu — 1, ^A = 0, Ho = 50 km s~^ 
Mpc~^. The luminosity is calculated in the observed 
frame ( no fc-correction) to compare with the ROSAT 
results (|Mivaii. Hasinger fc Schmidt 1 I2000D . We have 
assembled a total of 1004 AGN with Lq. 5-2.0 keV > 
10'*^ erg s~^, /x > 2 X 10~^^ erg cm~^ s~-^, and 



r' < 22.5 (Figure EJ. The ChaMP AGN boost the 
numbers at high redshift and lower luminosity due to 
Chandra' s faint limiting flux. Following the analysis of 
^iyaii. Hasinger & Schmidt (2000), we measure the co- 
moving space density for luminous (log Lq. 5-2.0 keV > 
44.5) AGN to ensure that we are sensitive to these ob- 
jects out to z ^ 4. This minimum luminosity ensures 
that our measurements are not highly biased by our X- 
ray flux limit. 

In Figure |51 we plot the co-moving space density of 
the 217 most X-ray luminous AGN. With these highly 
luminous AGN, we can compare the space density to past 
7? 05^1" results and optical surveys out to high redshift 
(z < 5). To easily compare the space density evolution 
from optical surveys, in Figure El we have renormalized 
their curves to match the ChaMP at z = 2.5. For all but 
the COMBO-17 survey, optical space densities {Mb < 
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Fig. 8. — X-ray luminosity vs. redshift of 1004 AGNs selected 
in the soft (0.5-2.0 keV) band with ROSAT and Chandra. The 
horizontal line shows our threshold for measuring the space density 
of highly luminous AGN. 
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Fig. 9. — Co-moving space density of 217 Chandra + ROSAT 
AGN selected in the soft (0.5-2.0 keV) band with log Lx > 44.5 
compared to the optical surveys. The optical space densities have 
been scaled to match the X-ray points at 2 = 2.5 for ease of com- 
parison. 



—26) have been scaled up by a factor of 3.2. COMBO- 
17 survey, which reaches much greater depths (R < 24) 
was scaled down by a factor of 7.1. Some caution must 
be taken when comparing to optical surveys with vastly 
different magnitude limits. For example, the AGN in the 
2dF survey evolve faster than those from X-ray selected 
samples at z < 2 (Figure 13. The COMBO-17 survey, 
which probes similar absolute magnitudes (Mb < —22.2) 
as the ChaMP, may show evolution rates comparable to 
the X-ray surveys over the redshift interval 1.5 < z < 2.1. 
The number density of our combined sample is similar 
to the ROSAT results with z < 3. Our larger z > 3 
sample (12), while statistically consistent with the ear- 
lier results of iMiva ii. Hasinger & Schmidt (2000), has 
smaller errors and shows a decrease in the space density 
from 2<z<3to3<z< 4.2 significant at the 4.6cr 
level. Even with the addition of ROSAT AGN, we have a 



slightly smaller sample at z > 3 than in the broad band 
sample (Section I6.1|l due to lower fluxes, luminosity or 
counts in the soft band for the ChaMP AGN. We find 
reasonable agreement with a slight excess between the 
relative drop in the space density of the X-ray and op- 
tical surveys at z > 2; our data are consistent with a 
decline in the X-ray luminosity function similar to that 
observed in the optical. 

7. CONCLUSION 

We have measured the broad band (0.3-8.0 keV) co- 
moving space density using a sample of 368 X-ray emit- 
ting AGN detected by Chandra. Our spectroscopic mag- 
nitude limit to date allows the inclusion of optically ob- 
scured AGN to z < 1 and unobscured AGN out to z '^ 4. 

Our primary results are as follows: 

• We confirm that low luminosity AGN (log Lx < 
44.5) are more prevalent at z < 1 than highe r 
luminosity AGN, as seen bv'Cowi e et al. I l)2003j) . 
[Fiore et al. (2003), and Ucda et al. I l|20fl3D . Non- 
BLAGN are the major contributor to the co- 
moving s pace density at these low redshifts, as re- 
ported bv lSteffan et al. I (|1)03). 

• The space density of type 1, highly luminous 
(log Lx > 44.5) X-ray selected AGN rises from 
the present epoch to a peak at z ~ 2.5 and then 
declines at z > 3. This behavior is similar to that 
of optically selected surveys. This is the first X-ray 
selected survey to detect a decline at high redshifts. 
The evolution is evident using a sample of 311 AGN 
from the ChaMP selected in the broad (0.3-8.0 
keV) Chandra band and 217 luminous QSOs se- 
lected in the soft band (0.5-2.0 keV) from Chandra 
and ROSAT. 

Our results support a more rapid depletion of fuel for 
the hi gh luminosity AGN from z ~ 2.5 to the present 
epoch l|Menci et al. 11200^ . perhaps due to their higher 
accretion rate. At z < 1, the lower luminosity AGN 
evolve slowly, which may be attributable to a low ac- 
cretion rate, compared to the more luminous AGN. At 
z > 3, the decline in the space density of highly luminous 
AGN represents the growth phase of supermassive black 
holes during a period of rapid galaxy assembly. The un- 
certainties in the X-ray luminosity function will improve 
with larger samples and a higher fraction of source clas- 
sification, either through spectroscopic or photometric 
techiques. With the inclusion of obscured, highly lumi- 
nous QSOs, X-ray surveys are well on the way to pre- 
senting a more comprehensive view of AGN evolution. 
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